Deterministic entanglement and tomography of ion spin qubits 
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We have implemented a universal quantum logic gate between qubits stored in the spin state of 
a pair of trapped '*''Ca ions. An initial product state was driven to a maximally entangled state 
deterministically, with 83% fidelity. We present a general approach to quantum state tomography 
which achieves good robustness to experimental noise and drift, and use it to measure the spin state 
of the ions. We find the entanglement of formation is 0.54. 



In recent years, entanglement has been detected and 
used in a variety of physical systems, but the number 
in which it can be created under deterministic control 
with good fidelity remains small. We present the imple- 
mentation of a universal 2-qubit quantum logic gate in 
a trapped ion experiment, demonstrating deterministic 
entanglement with 83% fidelity. Tomography was per- 
formed using a method robust against noise sources of 
practical relevance in the experiments. Our results illus- 
trate the relative insensitivity of the gate to the initial 
temperature of the ions. 

The logic gate exploits the oscillating force method 

1^ IJi IHl ' O^'' work differs from previous trapped 
ion experiments in the physical nature of the qubit, the 
trapping parameters, the light fields used to implement 
logic at the ions, and the laser sources, as follows. 

First, each qubit is stored in a pure spin-half system, 
the spin of a ^°Ca ion its ground state. Previous work 
used hyperfine structure ^ i jS| or a pair of electron 
orbitals 0, flil llll | . The spin- half system is convenient in 
its simplicity. The absence of 'spectator levels', i.e. other 
states nearby in energy, reduces the movement of popu- 
lation out of the controlled Hilbert space. This process, 
sometimes called leakage error, is not directly correctable 
by quantum error correction [l^ H^ . Spectator levels 
also exacerbate decoherence by photon scattering during 
the gate operation. 

Next, our trap is comparatively large Q, so electric 
field noise from fluctuating patch potentials and Johnson 
noise in the electrodes is very small 0|. We have 
measured very low heating rates, of order a few phonons 
per second, in the trap. 

We use a single laser field to apply all the operations 
(single ion rotations and the 2-qubit gate), changing only 
the amplitude and frequency to switch from one opera- 
tion to another. This reduces the experimental complex- 
ity, and opens the possibility of future work in which 
the different parts of a pulse sequence could be brought 
together in a single chirped and shaped pulse. 

Finally, the laser sources in our apparatus are all small 
semiconductor diodes, which do not require frequency 
doubling and which could in principle be packaged on a 
small optical chip. We thus begin to address the optical 



part of the challenge to scale this type of system up to a 
powerful quantum computer|l7l Il8lll9||. 

Let the spin-up/down states |T) j |i) be the computa- 
tional basis states of a qubit. We implement the 2-qubit 
controUcd-phase gate Zi((/)i)Z2((/)i)G'(^) where Zk{(j)i) is 
a rotation of qubit k — \,2hy (j)i about the quantization 
axis, and G'(^) is a two-qubit operator represented in the 
computational basis |tt) , ITi) i UT) : Ui) by the diagonal 
matrix diag(l, exp(i5'), exp(i\I'), 1). For = 7r/2, or an 
odd multiple of 7r/2, this is equivalent to controlled-not 
up to single qubit rotations. 

The gate mechanism is described in H. A spin- 
dependent force oscillating at frequency u is applied to 
the pair of ions. Let fkim) be the complex amplitude of 
the force on ion k when it is in spin state m ='\ , |. The 
COM (stretch) mode is excited by the sum (difference) 
force /i(toi) ± /2(to2) respectively. The effect on the 
quantum harmonic motion is simply to displace the state 
in its z~p phase space (as long as the position-dependence 
of the force is negligible during the displacement ^3)- 
By choosing uj close to one of the normal mode frequen- 
cies one can enhance the displacement and also simplify 
the dynamics, because then the excited mode simply de- 
scribes a circle in phase space, which closes after time 
2t:/S, and the excitation of the other mode can be ne- 
glected. This motion causes the system to acquire a 
phase (f>L proportional to the area of the loop in phase 
space. 

Many arrangements are possible. We adopt /i(m) — 
/2(to), then the stretch mode is excited only for anti- 
aligned spin states. By tuning close to the stretch mode 
frequency tOs, i.e. \d\ <C lOs — uJc where d = cu — LOg 
and LOc is the COM mode frequency, we then have 
^ = .Pl = (7r/2)(%/<5)2 where % = |/(T) - /(i)|zn 
zos = (ft/ (4Mci;s))^/^ and M is the mass of one ion 

To measure the spin state of the ion pairs, we devel- 
oped the following quantum state tomography method. 
Tomography in general requires the accumulation of a 
sufficient set {Mi} of measurements of a given state p 
(many copies of which are needed), in order to allow p 
to be reconstructed to some desired accuracy from {AU}. 
Usually rotations are applied to the system which is then 
measured in a fixed basis. 
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This problem has been studied in various settings and, 
for example, sets of {Mi} of minimal size, or designed 
for certain experimental situations, have been discovered 
[2^ . Our method is designed to be robust against noise 
and drift problems, and to be convenient in a fairly broad 
range of experimental settings. We achieve this by a judi- 
cious choice of rotation angles and by accumulating data 
sets in a form that should be well-fitted by sums of or- 
thogonal sinusoidal functions of known period. The least 
squares fit is robust because the functions in the fit are 
orthogonal, it averages the data over many timescales, 
and many noise processes affect the residuals not the fit- 
ted parameter values. 

Let R{9, 0) be a rotation of a single qubit through 9 
on the Bloch sphere about an axis of azimuthal angle 4>. 
To analyse a given density matrix p, the qubits are first 
rotated and then measured in the |t) , |i) basis, then p 
is re-prepared and one repeats with either the same or 
new rotation angles, until a sufficient set {Mi} has been 
acquired. Experimentally 9 depends on a pulse area, 
on the phase of an r.f. oscillator. The latter is more 
easily adjusted than the former to a variety of precisely 
known values. Therefore we only use two values of 9, 
but many values of 0, see fig. 2. For clarity we give 
details for the case of two qubits both undergoing the 
same rotation, but the principles apply more generally. 
In this restricted case a general p can be almost, but not 
fully reconstructed. 

Write p = J2'i j^o'^ij'^i ^ where Oi are the Pauli 
spin operators (starting with ao = I the identity). Let 
P{mim2) be the populations of p and Pj^imi, TO2) be the 
populations after the rotation R{9, (p) has been applied 
to both qubits. Then 

-PfidT) = a|| + fell cos((/)) + C|| sin((/)) 

-|-(i|| cos(20) + e|| sin(20) (1) 

and similarly for tiiiti where the coefficients a-e are 
simple sinusoidal functions of 9. By choosing different 
values of 9 and (p, 12 independent real numbers can be 
extracted from the population measurement outcomes. 
Therefore 3 of the 15 independent real numbers which 
fully characterize p are not available when both qubits 
undergo the same rotation. 

The tomography method consists of accumulating data 
at given as a function of (/), then fitting eq. ^ and its 
partners for Pr{\[)^ Pr{['\) to the data, with the a-e co- 
efficients as fitted parameters, c.f. fig. 2. After using two 
9 values, the density matrix coefficients Cij can be found 
from 9 and a-e (with some redundancy). The inferred 
matrix p^ is not guaranteed to be positive definite. One 
can handle this by any suitable approach, for example 
by searching among physical density matrices p^ for the 
one which most closely matches p*^ by some measure. 
We adopt this "maximum likelihood" method using the 



The apparatus consists of a linear r.f. Paul trap and 
laser system largely as described in We load single 
or pairs of ''^"Ca ions by photoionisation. A 6.3 MHz r.f. 
field provides radial confinement, with radial COM fre- 
quency ~ 1 MHz. Axial confinement is provided by d.c. 
endcaps, tuned as described below. 397 nm and 866 nm 
lasers are used for Doppler cooling, fluorescence detection 
and spin state preparation (with 99% fidelity) by optical 
pumping. A further pair of lasers performs spin-state- 
selective shelving to -D5/2, which allows the fluorescence 
to be used to read out the spin state (with '-^ 90% fi- 
delity), see 131. The ions are cooled in 3 dimensions by 
Doppler cooling, and along z by Raman sideband cooling 
to ncom.str — 0.2,0.2, measured as in ref. [isf . 

The light field used for pulsed sideband cooling and all 
the coherent operations is a walking wave created by laser 
beams A, B propagating at 9l = 58.9° to each other, 
with difference wavevector AA: = 47rsin(0j;,/2)/A directed 
along the trap axis z. A quantization axis is set by a weak 
magnetic field directed at 3° to beam B (and 9a = 62° 
to beam A). All these axes are horizontal. To reduce 
sensitivity to pulse area we balance the a.c. Stark shifts 
on It), I J,) from either beam acting alone, using Ramsey 
interferometry. This results in polarization very close to 
linear j2^ , which we make at angle (3 to the vertical for A, 
and horizontal for B. Each beam illuminates both ions 
equally. They are derived from the same laser, at a fre- 
quency detuned Ai/27r = 30 GHz above the S1/2 — P1/2 
transition. A precise frequency difference uj is introduced 
between them by acousto-optic modulation. 

The Zeeman splitting of |t) from ||) is (jJq = 27r x 
4800 kHz. At bj — ujQ spin-flip transitions are resonantly 
driven by the tt component of A combined with the a~ 
component of B. Let fJc be the Rabi frequency for this 
'carrier' Raman process. The oscillating force is created 
by the periodic light shift from the ct* components of 
the walking wave. The coupling strength fi^ is given by 
= ^/^^c^+l ^-m where e+,e7r are polarization ampli- 
tudes. For our geometry, \/2e+/e-,^ — coi9A/ cos(A0/2) 
and Vlf = 2r/f2s sin(A0/2), where 77 = Afczos is the 
stretch mode Lamb-Dicke parameter, and At/i is the 
phase angle between the forces, given by tan(A(^/2) — 
(cos 9a tan/3)~^. 

During the force pulse at w ~ 0;^, the carrier process is 
driven off- resonantly at detuning ± cj. For i ^ 3> 
ilc, the qubit states are pushed apart by a light shift 
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This is of order 4 kHz in our experiments, and it gives 
rise to the single-qubit rotations Zfc(0i) in the gate, with 

(pi = AcT. 

To achieve fi{m) — /2(™) we adjust ujc so that the 
ions' separation d ~ 9 /xm is an integer number p of 
standing wave periods, A/cd = 2'np. We used p = 
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22, ijjc/^n = 500 kHz [r] = 0.133) in one experiment, and 
p = 21, u;c/2tt = 536.5 kHz (77 = 0.128) in another. 

A given experimental sequence consists of cooling, spin 
preparation in then a spin echo sequence with the 
force pulse W of duration r in one or both of the gaps, 
followed by an analysis pulse and then spin state mea- 
surement; see inset to fig. lb. The analysis pulse is the 
rotation R{0, (j)) in the tomography scheme. A perfect 
implementation would produce the maximally entangled 
state 
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before the analysis pulse, with r = 2(pi — tt/2. To as- 
sess the state p obtained in practice we use the fidelity 
F = maxr {E{r)\ p\E{r)). This compares p with the 
most closely matching member of the class \E{r)) 27]. 
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FIG. 1: a) Measured populations Pfl(TT), Pfl(Ti) + Pfl(iT) (•, 
□ ) vs. r after a spin-echo sequence with VV(r) in the first gap, 
at S/2n = 12.6kHz. Each point is the average of 500 repeats of 
the experimental sequence. The lines are fitted curves, using 
a model which assumes each qubit decoheres independently at 
rate F. b) Parity signal vs. (p after a 6 — 0.467r analysis pulse, 
for r = 77 /is. Each point is the average of 1000 repeats of the 
experimental sequence. The inset shows the pulse sequence. 



We use "Schrodinger cat" experiments with a single 
ion [20I to analyse the forced motion and charac- 
terise light shifts. We then use the two-ion behaviour as 
a function of r with no analysis pulse to diagnose the 
setup. We model the results by assuming each qubit de- 
coheres independently at rate F. For the case of a W 



pulse only in the first gap we then expect 

P(tT) = A- ie-r"-l"(")l'/2cos(*(r))cos(A,T) 
Pill) + mt) = 1-2A (4) 

where A = (1/4) + e-^^''[cos{2AcT) + e-2|"(^)l']/8 and 
^(t), ^'(t) are the motional displacement and phase as 
described in ref Data is shown in fig la, fitted with 
floated parameters T,S,ilf and A(f). We obtained F — 
5.4ms-\ 6/2Tr = 13kHz, %/27r = 23kHz and = 1.6 
from the fit; these values were consistent with our other 
information on 6, fie and A0. Note that Qf/5~ V3 so 
in this experiment ^' = 37r/2. When we use the model 
and fitted values to infer p at r = 27r/(5, we obtain fidelity 
F = 0.7(1). 

We next measured a lower bound more directly, F > 
2\C\, by using an analysis pulse at 6 = tt/2 to deduce the 
coherence C = Pii,^ = cn — C22 + *(ci2 -|- C21). This is 
obtained from the component at frequency 2 in the parity 
signal PR(tJ,)+PR(|t) vs. (data shown in fig. lb.). We 
observed F > 0.74(3). This result is more precise than 
the previous one because it does not rely on assumptions 
about the decoherence. 
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FIG. 2: Tomography signals for spin state after the spin echo 
(W(44Ats) in both gaps, S/2-n- = 22.7kHz). The data points 
showPfl(TT),Pfl(n)+J'fl(iT) (•, □)vs. for 6» = 0.667r (top) 
and — 0.547r (bottom), as deduced from 500 repeats of the 
experimental sequence. The lines are the fitted curves Q. 
values were deduced from carrier flopping signals, taking into 
account a 0.1 /is dead time in the AOM; (p is accurately known 
from the r.f. signal generators. Inset: locus of 9, (j) values on 
the Bloch sphere. 



Finally, we performed tomography in another exper- 
iment where the W pulse was inserted in both gaps. 
This has the effect of cancelling the single-bit rotations 
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FIG. 3: Density matrix obtained from tomography data 
shown in fig. 2. The height of the bars indicates the abso- 
lute value of each density matrix element, and the clock faces 
indicate their phase. For this example the density matrix is 
closest to the state (3) with r — I.IStt. 



Zk{(j>i)- Both pulses had duration r = 2tt/6 so the mo- 
tion completes 2 loops in phase space, and the gate re- 
quires ^' — tt/A for each loop. We set 6/2tt ~ 22.7 kHz 
and observed the signal with pulses of double length in 
order to adjust the laser intensity. We then used the car- 
rier flopping rate to infer f2/, obtaining 27r x 16.3(9) kHz, 
giving the consistency check (f2//(5)^ — 0.52(5). 

The data and fitted curves for the tomography are 
shown in fig. |21 and the inferred (maximum likelihood) 
density matrix is shown in fig. O Owing to the absence 
of single ion addressing in the rotations and the measure- 
ments, the tomography only gives partial information on 
the internal elements of p (i.e. away from the corners in 
fig. OJ. However, we find that there is negligible pop- 
ulation in and ||t) and therefore for the state in 
question, the tomography is complete. The inferred p 
had fidelity F — 0.83(2) and entanglement of formation 
|l3| 0.54. The improvement is owing to better precision 
in pulse areas and a reduction in the total area by 
which reduces decoherence by photon scattering. 

We estimate that the primary source of infidelity is 
photon scattering (12%), with additional contributions 
from imbalance in light intensity at the ions, pulse area 
imprecision, and motional decoherence. Had we used a 
gate method which is more sensitive to the prepared mo- 
tional state, such as [l^, the thermal effect alone would 
contribute a further 6% infidelity through imprecise side- 
band pulses. We therefore obtained a useful improvement 
by using a temperature-insensitive method. 
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